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Abstract 


We present the final results of our search for new Milky Way (MW) satellites using the data from 
the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey over — 1,140 deg?. In 
addition to three candidates that we already reported, we have identified two new MW satellite 
candidates in the constellation of Sextans at a heliocentric distance of Do ~ 126 kpc, and Virgo 
at Do œ 151 kpc, named Sextans Il and Virgo III, respectively. Their luminosities (Sext Il: My ~ 
—3.9 mag; Vir III: My ~ —2.7 mag) and half-light radii (Sext II:rj, ~ 154 pc; Vir Ill:r;, ~ 44 pc) place 
them in the region of size-luminosity space of ultra-faint dwarf galaxies (UFDs). Including four 
previously known satellites, there are a total of nine satellites in the HSC-SSP footprint. This 
discovery rate of UFDs is much higher than that predicted from the recent models for the 
expected population of MW satellites in the framework of cold dark matter models, thereby 
suggesting that we encounter a too many satellites problem. Possible solutions to settle this 
tension are also discussed. 


© 2018. Astronomical Society of Japan. 
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1 Introduction 


The nature of dark matter in the Universe is still enigmatic 
(e.g., Bullock & Boylan-Kolchin 2017; Ferreira 2021). For 
many years, cold dark matter models in a A-dominated 
Universe (ACDM) have been the theoretical paradigm, 
because the models are successful for reproducing the 
observed large-scale structure of the Universe, at scales 
larger than ^ 1 Mpc (e.g., Tegmark et al. 2004; Planck 
Collaboration VI 2020). The theory is based on the hi- 
erarchical assembly of smaller to larger dark halos, be- 
cause the amplitudes of the initial density fluctuations are 
larger at smaller scales. Thus, smaller halos have higher 
internal densities and so were able to collapse earlier in 
the expansion of the Universe. Merging and accretion of 
these smaller halos in a hierarchical manner form more 
massive, larger halos as well as larger cosmic structures, 
which are found to be in agreement with various observa- 


tions of structures in the Universe. 


However, this assembly process leads to the survival 
of dense parts of the dark matter distribution, giving 
rise to hundreds of small subhalos surviving in a Milky 
Way (MW)-sized host halo, in apparent contradiction with 
the observed population of only tens of MW satellites 
(see, e.g., Bullock & Boylan-Kolchin 2017; Drlica- Wagner 
et al. 2020) — some of them are known as "classical" 
dwarfs with V-band absolute magnitude, Mv, in the range 
—18 mag € My < —8 mag, which were discovered before 
the Sloan Digital Sky Survey (SDSS) (York et al. 2000). 
This is called the missing satellites problem (Klypin et 
al. 1999; Moore et al. 1999), which has been a key chal- 
lenge to keep ACDM models. Also, denser smaller halos, 
which collapsed earlier, tend to sit in the central parts 
of a larger halo after merging takes place, producing a 
cuspy central density, whereas some or many of observed 
dwarf spiral and spheroidal galaxies show a core-like cen- 
tral density in their dark matter distribution (Oh et al. 
2011; Hayashi et al. 2020). This is called the core/cusp 
problem. Due to these small-scale issues in ACDM models, 
say at scales smaller than ~ 1 Mpc, several alternative dark 
matter models, which suppress the amplitude of the power 
Spectrum at smaller scales, have been proposed, such as 
Warm Dark Matter (WDM), Fuzzy Dark Matter (FDM), 
and Self-Interacting Dark Matter (SIDM) (e.g., Spergel & 
Steinhardt 2000; Bullock & Boylan-Kolchin 2017; Hui et al. 
2017; Ferreira 2021; Chan et al. 2021; Kondo et al. 2022). 


Among these small-scale issues in ACDM models, this 
work attempts to provide new clues to understand the 
missing satellites problem. Namely, there is still a pos- 
sibility that we are undercounting the true population of 
MW satellites, especially ultra-faint dwarfs (UFDs) with 
V-band absolute magnitude, My, fainter than —8 mag, 
due to various observational biases related to the survey 
area and depth (Koposov et al. 2008; Tollerud et al. 2008). 
However, this limitation in our knowledge of MW satel- 
lites has been overcome by recent massive imaging sur- 
veys, including the SDSS, the Dark Energy Survey (DES) 
(Abbott et al. 2016), and the Pan-STARRS 1 (PS1) sur- 
vey (Chambers et al. 2016). Indeed, these modern surveys 
have discovered a number of new UFDs, giving important 
insights into the missing satellites problem (e.g., Willman 
et al. 2005; Sakamoto & Hasegawa 2006; Belokurov et al. 
2006a; Laevens et al. 2014; Kim et al. 2015; Kim & Jerjen 
2015; Laevens et al. 2015a; Laevens et al. 2015b; Bechtol 
et al. 2015; Koposov et al. 2015; Drlica-Wagner et al. 
2015; Homma et al. 2019). These survey also imply that 
faint, diffuse UFDs at larger heliocentric distances, i.e., 
in the outer parts of the MW halo, still remain unde- 
tected and these are actually accessible with Subaru Hyper 
Suprime-Cam (HSC, Miyazaki et al. 2018). 

The purpose of this paper is to report the final results 
of our extensive search for UFDs from the HSC Subaru 
Strategic Program (HSC-SSP) (Aihara et al. 2018a; Aihara 
et al. 2018b). 
1.5 deg diameter field of view, and a 330-night survey with 
this camera has been conducted in the HSC-SSP frame- 
work. We have already reported the discovery of three 
new UFDs: Virgo I, Cetus III, and Boótes IV, from the 
internally released data before 2018 April (818A) over 
an area of ~ 676 deg? (Homma et al. 2016; Homma et 
al. 2018; Homma et al. 2019). Here, using the data ob- 
tained before 2021 January (internal data release 821A), 


This prime-focus camera, HSC, has a 


we present the discovery of two more candidate UFDs, 
Sextans II and Virgo III, as well as the identification of 
seven previously discovered or known satellites. In to- 
tal, there are nine satellites in the HSC footprint, i.e., five 
new satellites discovered by our group in addition to four 
previously known satellites: Sextans, Leo IV, Leo V, and 
Pegasus III. This discovery rate sets important constraints 
on the missing satellites problem in ACDM models. 

This paper is organized as follows. The data of HSC- 


SSP and the method to derive stellar overdensities using 
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The results of 
searching for new satellites in the HSC-SSP footprint are 


isochrone filters are given in Section 2. 


presented in Section 3. Section 4 is devoted to the discus- 
sion on the implications for dark matter models, and our 
conclusions are presented in Section 5. In this paper, mag- 
nitudes are given in the AB system (Oke & Gunn 1983). 


2 Data and Method 


The HSC-SSP survey provides multi-band (g, r, i, z, y) 
photometry data in the northern sky (Aihara et al. 2018a; 
Aihara et al. 2018b; Furusawa et al. 2018; Kawanomoto 
et al. 2018; Komiyama et al. 2018; Miyazaki et al. 2018). 
In the Wide layer that we adopt here, the sky coverage is 
~ 1,200 deg? and its target limiting magnitudes for a 5c 
point source are given as (g, r, i, z, y) — (26.5, 26.1, 25.9, 
25.1, 24.4) mag. As in our previous papers, we make use 
of the g, r, and i-band data obtained before 2021 January 
(internal data release S21A), which cover ~ 1,140 deg? 
in two separate fields along the celestial equator and one 
field around (a,6o) = (242?,43?) (Figure 1). These photo- 
metric data are processed using hscPipe v8.4 (Bosch et al. 
2018), which is a branch of the Legacy Survey of Space and 
Time (LSST) pipeline (see Ivezié et al. 2019; Axelrod et al. 
2010; Juric et al. 2017) calibrated against PS1 photometric 
and astrometric data (e.g., Tonry et al. 2012; Schlafly et 
al. 2012; Magnier et al. 2013). These HSC-SSP data are 
corrected for the mean foreground extinction in the Milky 
Way (Schlafly & Finkbeiner 2011). 


2.1 Selection of target stars 


'The method for the selection of target stars is the same as 
that adopted in our previous work: (1) point-like images 
are selected to avoid galaxies with extended images, (2) 
the criterion, g — r « 1.0, is placed to avoid a foreground 
contamination from M-type main-sequence stars, and (3) 
to further separate star candidates against other contami- 
nants, the fiducial g — r vs. r — i relation expected for stars 
is used. We describe each of these steps in more detail in 
the following. 

First, we select point sources based on the extendedness 
parameter calculated by the pipeline, which is defined in 
terms of the ratio of PSF to cmodel fluxes (Abazajian et al. 
2004), fesr / femodel! an object with fesr/ femodel > 0.985 
is regarded as a point source (extendedness= 0). Here, 
we adopt this parameter measured in the i-band, since the 
typical seeing in this band is the best of the five filters, with 
a median of ~ 0.6". Also, we decide to select point sources 
down to i = 24.5 mag as stars, because the completeness of 


stars in the star/galaxy separation is larger than 9096 at 


i < 22.5 mag, and decreases to about 50% at i = 24.5 mag, 
based on comparing the HSC data in the COSMOS field 
with that from HST/ACS (Aihara et al. 2018b). 

We then set the color criteria to avoid foreground disk 
stars and other contaminants, i.e., background quasars and 
distant galaxies, because these contaminants may often ap- 
pear as point sources. For this purpose, we utilize the field 
called WIDE12H at (RA, DEC) = (180°,0°), where the 
star/galaxy separation is robust for bright objects with 
i < 21 mag, and adopt the g—r vs. r—i color-color di- 
agram for stars with extendedness= 0 and galaxies with 
extendedness= 1. As shown in Homma et al. (2018) (Figure 
1), stars are distributed in a narrow sequence in this color- 
color diagram. We first remove M dwarfs in the MW disk 
characterized by the color g — r > 1, then set color cuts 
to select remaining candidate halo stars, with boudaries 
(g—r, r—i)=(1.00, 0.27), (1.00, 0.57), (—0.4, —0.55), and 
(—0.4, —0.25) and with a width of A(r — i) — 0.3 mag, 
which is wider than twice the typical 1o photometric error 
of r — i at i= 24.5 mag (see Homma et al. 2018). We utilize 
these color cuts to select the MW halo stars in this work. 


2.2 Algorithm for the detection of candidate UFDs 


As in our previous work (Homma et al. 2016; Homma et al. 
2018), we search for overdensities of stars compared to ran- 
domly distributed field stars and distant galaxies/quasars, 
after filtering stars using the isochrone filter as shown be- 
low (Koposov et al. 2008; Walsh et al. 2009). 


2.2.1 Isochrone filters 
Our isochrone filters for finding UFDs are designed to se- 
lect metal-poor, old stars inside a color-magnitude diagram 
(CMD), since the known UFDs in the MW are character- 
ized by the stellar populations as found in the MW glob- 
ular clusters. These filters are matched to finding such 
stars inside a CMD locus within a specified distance range 
from the Sun. The PARSEC isochrones (Bressan et al. 
2012) are adopted, with which we use four different mod- 
els, characterized as (a) an age of t= 13 Gyr and metal- 
licity of [M/H]= —2.2, (b) t = 13 Gyr and [M/H]= —1.5, 
(c) t=8 Gyr and [M/H]= —2.2, and (d) t 2 8 Gyr and 
[M/H]— —1.5 (See Fig. 2). 

These filters are defined in a CMD with g—r vs. r-band 
absolute magnitude, My. The width of each filter, as a 
function of r-band magnitude, is given from the quadrature 


sum of the lo error in the HSC photometry of the g— r 


color and a typical color dispersion (~ +0.05 mag) for red 


giant-branch stars (RGBs), which corresponds to a typical 


range of metallicities of ~ +0.05 dex within each of the 
observed MW UFDs. To search for UFDs in the distance 
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Fig. 1. Survey areas of the HSC-SSP Wide layer for the S21A release (bounded by red curves). The red and blue stars denote the locations of the newly 
found satellites in this work, Sextans II and Virgo III, respectively. The black squares and triangles denote, respectively, the known satellites and those found 
in the HSC-SSP reported earlier within these areas. We note that Pisces Il (shown with Pis2) is just outside the HSC-SSP footprint. 


18 


19 


20 


21 


22 


23 


24 


| 
jl \ | 


25 


EN Age:13Gyr, [M/H]=-2.2 
ENS Age:13Gyr, [M/H]=-1.5 
(a Age:8Gyr, [M/H]--2.2 
MN Age:8Gyr [M/H]=-1.5 


Fig. 2. Isochrone filters adopted in this work placed at distance modulus of (rn — M)o = 16.5, 20.0, and 24.0 mag, respectively. Four filters are shown with 
(1) an age of t = 13 Gyr and metallicity of [M/H]= —2.2 (green area), (2) t = 13 Gyr and [M/H]= — 1.5 (blue shading), (3) ¢ = 8 Gyr and [M/H]= —2.2 
(yellow shading), and (4) t = 8 Gyr and [M/H]— — 1.5 (red shading). The width of each shaded region corresponds to the typical photometric error in color at 


each r magnitude with an additional 0.05 mag added in quadrature.. 


range of Do = 20 to 631 kpc, we shift the isochrone filter 
in a CMD over the distance moduli of (m — M)o = 16.5 
to 24.0 mag in steps of 0.5 mag. Figure 2 shows our four 
isochrone filters placed at (rn — M)o = 16.5, 20.0, and 24.0 


mag, respectively. 


2.2.2 Method to identify overdensities as UFD 
candidates 

Armed with these isochrone filters, we search for statis- 

tically significant overdensities as UFD candidates in the 

following manner. 

First, we select and count the stars in 0.05? x 0.05? 
boxes in RA and DEC with an overlap of 0.025? in each 
direction. This grid interval, 0.05°, in this search is com- 
parable to a half-light diameter of ~ 80 pc for a typical 
UFD located at a distance of Do = 90 kpc, so that the 
UFDs at Do > 90 kpc, i.e., our targets in HSC-SSP, can 
be found within these grids. 


Second, we assemble the number of stars in each cell, 
ni, j (i for RA, j for DEC), where the cells with no stars, 
nij = 0, e.g., due to masking near a bright-star image, 
are ignored, and then estimate, for each of the separate 
Wide-layer fields, the mean density (7) and its dispersion 
(c) over all cells. Based on these values, we define the 
normalized signal in each cell, 5;,;, which is the number of 
standard deviations relative to the local mean (Koposov et 
al. 2008; Walsh et al. 2009), Si; = (ni, — n)/c. We find 
that this distribution of S is nearly Gaussian, as we found 
in our previous papers. 

Finally, to select overdensities as UFD candidates that 
have statistically high significance, so that false detections, 
such as random fluctuations, are avoided, we adopt a de- 
tection threshold, Sin, for the value of S (Walsh et al. 
2009). To obtain St, we performed a Monte Carlo analy- 
sis as reported in our previous paper (Homma et al. 2018), 


where for purely random fluctuations in stellar densities, 
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the maximum density contrast, Smax, as a function of n is 
derived. Then, we adopt a conservative choice for a thresh- 
old in the density contrast of Sth — 1.5 x Smax, with which 


we identify statistically significant overdensities as UFDs. 


3 Results 


We run the search algorithm explained above on the 821A 
data of the HSC-SSP. As a result, more than 30 over- 
densities exceeding the threshold are detected, but about 
half of them are found to be false due to contamination 
from artifacts, such as those around bright stars. We dis- 
tinguish between fake and real overdensities of stars, by 
looking at the gri color HSC image and the CMD around 
its location for each of the candidates. The other half of 
the overdensities are real stellar systems, such as globu- 
lar clusters(GCs): Pal 3, Whiting 1 (Harris 1996), star 
cluster(SC): HSC 1 (Homma et al. 2019), and dwarf ir- 
regular galaxies(dIrrs): IC1613, Sextans B, and DDO 190 
(McConnachie 2012). The significance and statistical val- 
ues defined in Section 2.2.2 of all the overdensities identi- 
fied as real stellar systems are given in Table 1. Note that 
the dIrrs are located outside the range of (m — M)o < 24.0 
and not associated with MW, although they are detected 
by their TRGBs (Tip of the red-giant branch stars) passing 
through the bottom of the ishochrone filter at (m — M)o ~ 
16.5. 

In this section, we summarize the nine MW satellites 
or candidates detected in the HSC-SSP footprint; all of 
these will be simply referred to as 'MW satellites’ in what 
follows unless otherwise noted. While four of these satel- 
lites, Sextans, Leo IV, Leo V, and Pegasus III, are al- 
ready known, the other five, Virgo I, Cetus III, Bodtes 
IV, Sextans II, and Virgo III, are candidate satellites iden- 
tified in the HSC-SSP data, among which Sextans II! and 
Virgo III are newly found in the S21A data and reported 
in this paper for the first time. 


3.1 Summary of the MW satellites in the HSC-SSP 
footprint 


We describe below the procedure to calculate the basic 
properties of these MW satellites, including their heliocen- 
tric distances, half-light radii, and absolute magnitudes. 
These quantities are summarized in Table 2 for all nine 
satellites identified in the HSC-SSP footprint, where we 
include the results for four known satellites, Sextans, Leo 
IV, Leo V, and Peg III, to confirm that our current analy- 


sis indeed yields comparable results obtained by previous 


! Parallel to our study and independently, Gatto et al. (2023) also reported 
the discovery of Sextans II in the Kilo-Degree Survey (de Jong et al. 2013). 


works. Three candidate satellites that we reported ear- 
lier: Vir I, Cet III, and Boo IV, are re-analysed for con- 
sistency with the results of our previous analysis (Homma 
et al. 2018; Homma et al. 2019); with the exception of 
the estimate of V-band absolute magnitudes (see below), 
all measured properties remain unchanged for these three 
satellites. 

Here, we identify two new candidate MW satellites, 
which we name Sextans II (Sext II) and Virgo III (Vir IIT) 
at (RA, DEC)= (156.4°, —0.6°) and (186.3°,4.4°), respec- 
tively. In Figures 3 and 4, we show the space and CMD 
distribution of stars in these candidates, demonstrating 
that they are true stellar systems: RGBs, main sequence 
turn-off stars (MSTOs), and blue horizontal-branch stars 
(BHBs) are clearly visible in the CMD. Vir III is found to 
have the smaller number of each of RGBs, MSTOs, and 
BHBs than does Sext II. These are considered to be UFDs 
because (1) their distributions in CMDs are qualitatively 
consistent with old stellar populations and (2) they are spa- 
tially extended compared to globular clusters with similar 
brightness, as we describe in Section 4.1. 

For completeness, in Appendix 1, we also show the 
space and CMD distributions of stars in the other seven 
satellites (Figures 10-16). We note that the methods to de- 
rive the structural parameters (Section 3.3) and absolute 
magnitude (Section 3.4) of these MW satellites are suit- 
able only for very low luminosity systems such as UFDs, 
in which the number of detected stars is as small as tens 
to hundreds. Thus, Sextans, a classical dSph with many 
more stars detected, is an exception for this analysis, and 


we only calculate its distance. 
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Table 1. Significance of the overdensities identified as stellar systems 


Name Classification Filter^ (m — M)o (à, j) Nig n c Sg Smar Si j/Smas 
[mag] [deg] 
Sextans dSph(Classical) d 20.0 153.250,—1.625) 157 1.53 0.80 195.7 7.8 25.2 
Leo IV dSph(UFD) a 21.0 173.250, —0.550) 31 1.84 1.01 28.9 7.1 4.1 
Leo V dSph(UFD) a 21.5 (172.775,2.225) 38 1.86 1.02 35.5 7.0 5.1 
Pegasus III dSph(UFD) a 21.5 (336.100, 5.425) 25 2.08 1.13 20.3 6.7 3.0 
*Virgo I dSph(UFD) a 20.0 180.025, —0.675) 16 1.84 1.01 14.0 7.0 2.0 
* Cetus III dSph(UFD) a 22.0 (31.327, —4.285) 13 1.92 1.07 10.4 6.9 1.5 
*Bootes IV dSph(UFD) a 21.5 (233.700, 43.725) 26 1.49 0.75 32.9 7.9 4.2 
*Sextans II dSph(UFD) b 20.0 156.425, — 0.600) 21 1.68 0.85 22.7 7.5 3.0 
* Virgo III dSph(UFD) b 20.5 (186.350, 4.450) 25 2.08 1.15 19.9 6.6 3.0 
Pal3 GC b 19.5 (151.375,0.075) 322 1.56 0.80 399.5 7.7 52.0 
Whiting 1 GC b 16.5 (30.750, —3.250) 139 2.96 1.68 80.8 6.0 13.5 
*HSC 1 SC c 18.0 (334.325, 3.475) 31 1.53 0.79 37.5 7.8 4.8 
IC1613 dirr b 24.0 (16.275,2.025) 384 1.48 0.75 509.1 7.9 64.4 
Sextans B dirr d 16.5 (149.950, 5.300) 53 1.64 0.87 59.3 7.5 7.9 
DDO 190 dirr d 17.0 (216.200, 44.525) 10 1.80 0.58 15.1 8.5 1.8 


^ Models of the isochrone filters described in Section 2.2.1. 


* Objects discovered in the HSC-SSP survey. Their classifications are not yet uncertain. 
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Fig. 3. Sextans II. (a) The stellar overdensity passing the isochrone filter of t = 13 Gyr and [M/H]= —2.2 at (m — M)o = 20.5 for the point sources satisfying 
i < 24.5, g — r < 1.0, and the color-color cuts expected for stars in the g — r vs. r — i diagram (see subsection 2.1). (b) The spatial distribution of the galaxies 
which pass the same isochrone filter and constraints as for the stars. An overdensity at the center of this plot is seen, possibly because of misclassification of 
faint stars within a cluster as galaxies in hscPipe. (c) The spatial distribution of the objects classified as stars around the overdensity. The green circles and 
black dots denote, respectively, the stars inside and outside the isochrone filter at the distance modulus of (m — M)o = 20.5. A solid red curve shows an 
ellipse with a major axis of r = 2.0r;, (rj, = 4.2’) and an ellipticity of 0.43, whereas dotted red circles show annuli with radii r = 10.0r;, and r = 10.2r;, from 
the center of the overdensity. (d) The CMD of the stars in the g — r vs. r relation, which are located within the solid red ellipse in panel (c). Blue circles denote 
the BHB stars which we will use for an additional distance estimate as described in the text. (e) The same as (d) but for field stars at 10.0r;, < r < 10.2r;, 
having the same solid angle. None of the features of a stellar population (RGBs, MSTOs, BHBs) is apparent in this plot. 
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Fig. 4. Virgo Ill. Same as Figure 3 except that the isochrone filter is at (m — M)o = 20.9. (c) The solid red curve shows an ellipse with a major axis of 
r = 2.0ra (rj = 1.0’) and an ellipticity of 0.29. 


Table 2. Properties of previously known and candidate MW satellites in the HSC-SSP footprint 


Name? (R.A., Dec)? (1,b)° (m—M)o% | De of D Nh Th) Mos, v) Ay“ 
[deg] [deg] [mag] [kpc] [deg] ('] or [pc] [mag] [mag] 


Previously known MW satellites 


Sextans (153.262, —1.615)  (243.498,42.272) — 19.8*02 9178 = = — = — 0.168 
Leo IV (173.243, —0.536) (265.452,56.515) — 21.0'02 158116 110:53  0.17t998 10310 2.6t%2 or 119575. —4.50t0:35 0.069 
Leo V (172.784,2.222)  (261.851,58.535) 21.3192  182Ł18 106113 0.372042 50t Tote or 62119 —3.96+937 0.073 

: F3 or 88t22 Ali 93. 0.404 


Pegasus III — (336.100,5.413) 


M 
N 
= 
[on 
IT 
oo 
ioi 


>= 


69.851, —41.822 20e. SSP: Quam. 4E quus 


Candidate MW satellites 


Virgo I (180.038, —0.681) (276.942, 59.578) 19.8*02 91+3 6278, 0.594912  18*3  1.8t93 or 47413  —0.90*'072 0.066 
Cetus III (31.331,—4.270)  (163.810,—61.133) 22.0702 251438 10178  0.76*006 1672 1.2193 0r90*322  —3.45*015 0.066 
Bootes IV — (233.689,43.726) (70.682, 53.305) 21.0102 20082 3*4 0.6419:03 124110 ety or 462138 vetat 0.067 
Sextans II (156.437, —0.631) (245.326,45.322) 20.502 «126417  —1678 0.431505  93*10 42105 or 154435 —3-91t056 0.177 
Virgo III (186.348, 4.441) (286.476,66.477) 20.9703. 151413 adr. 029701 25t] . 10'020r44"1$  —-2.69*029 0.054 


? Satellite's name 

P Equatorial coordinates (J2000) 
* Galactic coordinates 

d Distance modulus 


(m— M) 
© Heliocentric distance calculated as 10-5 ~~? kpc 


* Position angle from north to east 

8 Ellipticity defined as € = 1 — (b/a) with a the scale-length of the system along its major axis and b that along its minor axis. 
n Number of detected stars estimated to be members of the stellar system 

i Half light radius 

J V-band absolute magnitude 


* Integrated magnitudes are corrected for the mean Galactic foreground extinction, Ay (Schlafly & Finkbeiner 2011). 
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Table 3. Distance moduli derived from BHBs 


Name Number of BHB stars? (m — M)o* 
[mag] 

Sextans 195 19.77 + 0.20 (19.8 + 0.20 
Leo IV 8 21.00 + 0.08 (21.0 + 0.20 
Leo V 8 21.29 + 0.05 (21.3 + 0.20 
Pegasus III 5 21.53 + 0.06 (21.5 + 0.20 
Bootes IV 10 21.54 0.19 (21.6 + 0.20 
Sextans II 9 20.47 + 0.22 (20.5 + 0.20 
Virgo III 3 20.91 + 0.04 (20.9 + 0.20 


"Within 5 times the half light radius of the centroid. 
>The values in parentheses are those derived from isochrone 


shifting, and listed for comparison. 


3.2 Distance estimate 


The first step after identifying the candidate stellar sys- 
tem is to estimate its heliocentric distance, or distance 
modulus (m — M)o, by fitting the isochrone filter to the 
overdensity having the highest significance (i.e., S defined 
in Section 2.2.2), where (m — M)o is shifted in steps of 0.1 
mag around the value derived in detecting it. This evalu- 
ation is mostly consistent with another distance estimate 
below, and we assume that the systematic uncertainty in 
It is found that for four 
known satellites, Sextans ((m — M)o =19.8+0.2 mag from 
this work), Leo IV (21.040.2 mag), Leo V (21.3+0.2 mag), 
and Peg III (21.5+0.2 mag), our distance estimates are 


distance modulus is 0.2 mag. 


within 1c of previous estimates, (m — M)o = 19.67 + 0.10, 
20.94 + 0.09, 21.25 + 0.12, and 21.56 + 0.20 mag, respec- 
tively (McConnachie 2012 for the first three satellites, and 
Kim et al. 2015 for Peg III). 

The corresponding heliocentric distances to the newly 
found satellites, Sext II and Vir III, are 126*1? kpc and 
181715 kpc, respectively. 


An alternative estimate for the heliocentric distance 
comes from the BHBs, given their absolute magnitude cal- 
ibrated by Deason et al. (2011). The method to select 
BHBs from the data uses three colors, g — r, r — ài, and 
i— z, as described in Fukushima et al. (2019). We apply 
this procedure to the satellites which contain BHBs (ex- 
cept Virgo I and Cetus III), and Table 3 shows the results. 
'These distances from BHBs are found to be in agreement 
with those shown in Table 2 with a difference of a few 


percents. 


3.8 Structural parameters 


The next step is to calculate the structural properties 
of each stellar system following Martin et al. (2008) 
and Martin et al. (2016). 
(ao, 60,9,€,7n, Nx); (ao, ôo) for the celestial coordinates of 


There are six parameters 


the centroid of the overdensity, 0 for its position angle from 


(ens 
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Fig. 5. The density distribution of the stars in Sext Il and Vir IIl passing the 
isochrone filter given in Figure 3 and 4, in elliptical annuli as a function of 
mean radius. The uncertainties are based on Poisson statistics. The red 
line denotes the best-fit exponential profile within 5 ~ 6r» plus a constant 
representing the background. 


north to east, e for the ellipticity, ra for the half-light radius 
measured along the major axis, and N, for the number of 
stars within the isochrone and brighter than our magnitude 
limit, which belong to the overdensity. We then obtain the 
best-fit parameter set and their uncertainties by perform- 
ing a maximum likelihood analysis following Martin et al. 
(2008) for the stars within a circle of 5 — 6 times the in- 
ferred half-light radius. 

We find that for the known UFDs, Leo IV, Leo V and 
Peg III, for which the above method is applicable, the de- 
rived structural parameters are in good agreement with the 
previous results: Leo IV: (rp (pc),e) = (119*25,0.17* 0 08) 
compared to (127.8, 0.05) by Sand et al. (2010), Leo V: 
(62*19,0.37* 019) compared to (65133,0.52) by Sand et 
al. (2012), and Peg III: (88*25,0.37*011) compared to 
(79*35,0.46*0:13) by Kim et al. (2015). 

Applying this method to Sext II and Vir III, we find 
(ra (pc), £) = (154*35,0.43*007) and (44*12.,0.29*0-15), re- 
spectively. 

In Figure 5, we show the radial distribution of the stars 
passing the isochrone filter depicted in panel (d) of each 
of Figures 3 and 4, which is obtained from the count of 
the average density within elliptical annuli. The red line 
in the figure denotes the best-fit exponential profile rep- 
resented by the best-fit parameter set above. The radial 
distributions of other satellites are shown in Appendix 1 
(Figure 17). 


3.4 V-band absolute magnitude 


Finally, we estimate the absolute magnitude of each stellar 
system in a similar manner to that described in Martin et 
al. (2008). In summary, we extract stars randomly one by 


one and sum up their V-band luminosities from the stel- 


10 Publications of the Astronomical Society of Japan, (2018), Vol. 00, No. 0 


lar population model (Bressan et al. 2012) with an age of 
13 Gyr and metallicity of [M/H]— —2.2 until the number of 
stars brighter than the threshold reaches N., and it is our 
magnitude limit (i — 24.5 mag) converted to the absolute 
magnitude V using the formula given in Jordi et al. (2006) 
and the distance modulus (m — M)o. In this way, we per- 
form a Monte Carlo procedure of 10,000 cases and estimate 
the most likely value of My and its uncertainty including 
that of N, and (m— M)o. Here we choose three initial mass 
functions (IMFs) for the stellar population: from Kroupa 
(2002), Salpeter (1995), and Chabrier (2001). Unlike total 
stellar mass, the absolute magnitudes of stellar systems are 
quite insensitive to the IMFs, because the number of low 
mass stars has little effect on the total luminosity (Table 
4). We note that these Mtot,v for Vir I, Cet III, and Boo IV 
are somewhat brighter than what we derived in our previ- 
ous work (Homma et al. 2018; Homma et al. 2019), where 
the V-band absolute magnitudes for stars brighter than 


our magnitude limit (i — 24.5 mag) were presented. 


With this method, for Sext II and Vir III, we obtain 
the V-band absolute magnitudes of —3.91*027 mag and 
—2.60* 0-25 mag, respectively. 


4 Discussion 


4.1 Comparison with globular clusters and other MW 
satellites 


In order to determine whether both Sext II and Vir III are 
in fact MW UFDs, their structural properties are com- 
pared with those of MW globular clusters and known 
dwarf satellites. Globular clusters are known to be self- 
gravitating stellar systems having compact sizes and round 
shapes, whereas dwarf satellites are characterized by a dif- 
fuse, extended distribution of stars having an ellipitical 
shape. Also, in contrast to globular clusters, the line-of- 
sight velocity dispersions of members stars in dwarf satel- 
lites suggest that they are surrounded by a massive dark 
halo. 


Figure 6 shows the relation between the half-light radii, 
Th, and V-band absolute magnitudes, My, for MW dwarf 
satellites (squares) (McConnachie 2012; Laevens et al. 
2014; Bechtol et al. 2015; Koposov et al. 2015; Drlica- 
Wagner et al. 2015; Kim et al. 2015; Kim & Jerjen 
2015; Laevens et al. 2015a; Laevens et al. 2015b; Torrealba 
et al. 2016; Torrealba et al. 2018; Cerny et al. 2023) and 
globular clusters (dots) (Harris 1996). Sext II and Vir III 
are depicted as the red and blue stars with error bars, re- 
spectively. This figure indicates that these stellar systems 
have systematically larger ra than those of MW globular 


clusters having similar My, and that these are located in- 
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Fig. 6. The relation between My (mag) and ra (pc) for stellar systems. 
Dots denote globular clusters in the MW taken from Harris (1996). Filled 
squares denote the MW dSphs taken from McConnachie (2012), the recent 
DES work for new ultra-faint MW dSphs (Bechtol et al. 2015; Koposov et al. 
2015; Drlica-Wagner et al. 2015), and other recent discoveries (Laevens et 
al. 2014; Kim et al. 2015; Kim & Jerjen 2015; Laevens et al. 2015a; Laevens 
et al. 2015b; Cerny et al. 2023), whereas open triangles denote those found 
from the previous data release of the HSC-SSP (Homma et al. 2016; Homma 
et al. 2018; Homma et al. 2019). The red and blue stars with error bars cor- 
respond to Sextans II and Virgo III, respectively. The dotted lines are lines of 
the constant surface brightness, py = 26, 31, and 32 mag arcsec~?. 


side the locus of MW dwarf satellites. Also, their shapes 
are quite flattened (e = 0.3 ~ 0.4), like other dwarf satel- 
lites, and much more so than globular clusters (Harris 
1996; Marchi-Lasch et al. 2019). Thus, we conclude that 
Sext II and Vir III are new UFDs in the MW. 


Next, we compare the heliocentric distances of Sext II 
and Vir III with the distance distributions of MW globu- 
lar clusters and dwarf satellites. Figure 7 shows the rela- 
tion between My and Do for these stellar systems, where 
red and blue stars with error bars correspond to Sext II 
and Vir III, respectively. In this diagram, red and blue 
lines show the detection limits of SDSS and HSC, respec- 
tively. The former limit is taken from the SDSS complete- 
ness distance given as AR = 10672" Mv b") Mpc with 
(a*,b*) = (0.205,1.72) (Koposov et al. 2008). The latter 
limit for HSC, Hx is evaluated from HI m = 
10°-2(r,HSC—™r,SDSS) where Mr spss and m, nsc are the 
limiting point-source magnitudes in r-band for these sur- 
veys, given as Mm, spss = 22.2 mag and m, nsc = 24.7 mag, 
respectively (Tollerud et al. 2008). 
Mr usc is taken from our cut-off i-band magnitude of 


The latter value of 


i = 24.5 mag and typical r — i colors of stars in MW satel- 
lites of r — i = 0.2 mag (Homma et al. 2019). As is clear 
from this figure, both Sext II and Vir III are located at he- 
liocentric distances beyond the radial distribution of typ- 
ical globular clusters and are also beyond the detection 
limit of SDSS. 


Publications of the Astronomical Society of Japan, (2018), Vol. 00, No. 0 11 


Table 4. Estimation of M,.¢,y and M. from various IMFs 


Name Mtot,v Mx 
[mag] [Mo] 
Kroupa Salpeter Chabrier Kroupa Salpeter Chabrier 

Leo IV —4.50+0:35  —4.52+0:32  _449+03T | 74t? x103 14752 X10!  T8'72x 10° 

Leo V -3.964935 —4.05+9:38 3.974935 | 4.7tŁ3 x103 9172 «09 — 45 rox TOP 

Pegasus III | —4.11*022 — —408'03,  —407*02. | 5.1214x10? 9.920 x10? — 4.8-E L4 x 10? 

Virgo I -0.0- 972 0.83937 — — 9.847023 | Si T9 x10? GT S x10? gati x10? 

Cetus III -3.454949 -3.43104 83805 | 26; x dU  5.3+25x103 SATUS 10" 

Bootes IV | —5.34ł073 -5.381628  —530'025 | 16'0:x10*  20*05x10*. 1.6193 x104 

SextansII | —3.91*0 37. —3.98*0 30 — —3.88*0 35 | 43*11 x103  84*22 x10? 3.9778 x103 

Virgo III 269 ne -2.81404  —AT7T8 aes | Lote xl? “Sot?  1.4+07 x103 
-14 uM liocentric distance, and accretion time to those of LMC, to 
as [| A tse candidate " . . f include the impact of the LMC in the analysis. They then 
anl mr . . ; ] consider a detailed model of the galaxy-halo connection, 
L Conti | to associate satellite galaxies with subhalos in the simula- 
* " [ tions (Nadler et al. 2019), using abundance matching, and 
E $ [ taking into account satellite sizes, satellite disruption due 
* at to baryonic effects, galaxy formation efficiency, and orphan 
2L satellites. This model is then marginalized over these as- 
oL trophysical uncertainties in the fit, to the observed satellite 
2 populations by DES and PS1 surveys over the very large 


De [kpc] 


Fig. 7. The relation between My (mag) and heliocentric distance (kpc) for 
the systems shown in Figure 6. The red and blue lines indicate the detection 
limits of SDSS and HSC, respectively. 


4.2 Implication for A CDM models 


In this S21A data release of the HSC-SSP survey, two new 
satellites, Sext II and Vir III, have been discovered, in ad- 
dition to the re-identification of three satellites, Vir I, Cet 
III and Boo IV, which were reported earlier using the previ- 
ously released HSC-SSP data (Homma et al. 2016; Homma 
et al. 2018; Homma et al. 2019). As we have presented 
in this work, there are four more known satellites in the 
HSC footprint, Sextans, Leo IV, Leo V and Peg III, in- 
dicating the presence of nine satellites in total over the 
area of — 1,200 deg?. The implication of this discovery 
rate of satellites for ACDM models is discussed here, by 
adopting the latest studies of the expected population of 
satellites in ACDM models by Nadler et al. (2020) (see also 
Nadler et al. 2019; Drlica- Wagner et al. 2020; Dooley et al. 
2017; Newton et al. 2018). 

To obtain the underlying MW satellite population, 
Nadler et al. (2020) adopt a high-resolution cosmological 
zoom-in simulation for MW-mass host halos selected from 
the work of Mao et al. (2015). In particular, they select 
two MW-like halos, which are required to have an LMC 


analog companion galaxy, i.e., having a similar mass, he- 


area of ~ 15,000 deg”. Based on this model, they arrive at 
a total number of Mot = N(Mv < 0) = 220 + 50 satellites 
with ra > 10 pc within the virial radius of the MW. In 
what follows, we apply this model to our results using the 
HSC-SSP Wide layer. 

First, we make a correction for N( Mv) associated with 
the survey area of HSC-SSP over 1,141 deg?, under the 


assumption of an isotropic distribution of satellites. This 


area corresponds to a sky fraction of fo,usc — 0.028. Thus, 
the number of predicted satellites over this sky area is 
fo, usc. N (My). 

Second, we make a completeness correction associ- 
ated with the detection limit of the HSC-SSP survey, 
which can be expressed in terms of the completeness dis- 
tance, HI (Mv), beyond which a satellite with partic- 
ular My cannot be detected (Figure 7). We then eval- 
uate the completeness correction factor, frusc(Mv) = 
N(« RES, (Mv ))/Neot, following the work of Nadler et al. 
(2020) to obtain the expected radial distribution of satel- 
lites using their code available at GitHub?. The result is 
shown in Figure 8 (blue bands), from which we evaluate 
N(« r). The total number of satellites expected for the 
HSC-SSP data is then given as fr nsc fo,usc N(Mv ). 

Figure 9 shows the 
(N(My), black bands), 


counts for all satellites 


for the sky-coverage correc- 


2 github.com/eonadler/subhalo. satellite. connection. We thank 
Ethan Nadler for allowing us to use this code for this paper and for his 
advice in reproducing the results. 
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Fig. 8. Expected radial distribution of MW satellites reproduced from the 
work of Nadler et al. (2020) (blue bands), compared to the observed satellite 
distribution from HSC (green histogram). The red curve shows the Einasto 
profile fit to the radial distribution of satellites predicted by Newton et al. 
(2018). 
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Fig. 9. Expected number of MW satellites obtained from the work of 
Nadler et al. (2020) (N (My), black bands), sky-coverage corrected count 
(£fo,usc N (Mv ), blue bands), and corrected further for the magnitude limit 
(f£-, usc fo,Hsc N(Mv), red bands), where dark and light bands corre- 
spond to 68% and 95% confidence intervals, respectively. Green histogram 
shows the observed luminosity function obtained from the current work with 
HSC except Sextans(My = —9.3 mag, McConnachie 2012). 


tion (fo,nsc.N(Mv), blue bands), and for the combined 
correction (fy, usc fo;nsc.N(Mv), red bands), where dark 
and light bands correspond to 6896 and 9596 confidence 
intervals, respectively. For My < 0 (and ra > 10 pc, 
which is fulfilled for the observed satellites), the total 
number of satellites expected for the HSC-SSP footprint, 
fruscfo,uscN(My), is estimated as 3.9 0.9. This 


expected number is significantly smaller than the nine 


satellites actually observed, whose cumulative luminosity 
distribution is shown as the green solid line in Figure 9. 
Thus, instead of a missing satellites problem, this is a 
problem of too many satellites in the HSC-SSP footprint. 
We had come to a similar conclusion with the S18A data 
release (Homma et al. 2019). 


This tension is not solved even if we adopt a different 


model for the radial distribution of satellites. For instance, 


we adopt the Einasto profile, 
InN (r)/N_2 = —(2/o)(r/r-a)" — 1] , (1) 


as shown in the model of Newton et al. (2018), where « 
is a shape parameter, r_2 is the radius at which the log- 
arithmic slope is —2, and N_2 is the density at r = r. 5. 
'The concentration of this density profile is parameterized 
by c200 = R200/r—2, where Hoo is the virial radius. If we 
adopt o = 0.24, c200 = 4.0, and Hoo = 200 kpc (derived 
for a 1.0 x 10? Mc MW-mass halo (Newton et al. 2018)), 
we obtain the expected number of satellites in the HSC 


footprint as 4.3 + 1.0, which again falls short of the nine 
observed satellites. 

While this high detection probability of satellites in the 
HSC-SSP appears to solve the missing satellites problem 
in ACDM models, the remained tension with the predicted 
models for the population of MW satellites needs to be 
settled. A possible key is the observed radial distribu- 
tion of nine satellites in the HSC footprint as shown in 
the green line in Figure 8. Compared with the model re- 
sult (blue bands), the observed radial distribution is more 
extended, although we are limited by small-number statis- 
tics. However, taking into account that the Subaru/HSC 
deep imaging survey is more sensitive to faint, distant MW 
satellites than other existing surveys, it is possible that 
the true radial distribution of MW satellites is more ex- 
tended than that derived from the DES and PS1 surveys, 
so that the matching of the simulated subhalos with the 
observed MW satellites needs to be reconsidered. Indeed, 
in this more extended case, the predicted number of satel- 
lites could be larger than the case of the currently adopted 
model by Nadler et al. (2020), settling the tension with 
HSC-SSP survey counts (Kim et al. 2018). Also, related to 
this point, more numerical simulation data for MW-sized 
dark halos having an LMC analog are needed for getting 
the higher statistical significance in these results. 

Secondly, it should also be borne in mind that the de- 
tectability of satellites in the HSC-SSP is simplified by as- 
suming it to be a step function which changes from 10096 
to 096 beyond a particular radius calculated as the de- 
tection limit distance (Rcomp) at each magnitude. Note 
that the detection limit from Koposov et al. (2008), which 
we adopt for our completeness boundary, indicates the 
5096 detection rate only to objects with surface bright- 
ness uy < 30 mag arcsec”? and ignores the complicated 
issue that the detectability varies with a location in the 
surveyed area due to the crowdness of the sky (Doliva- 
Dolinsky et al. 2022; Qu et al. 2023). Given this fact, 
in reality there should be a greater dispersion in the ex- 
pected number of MW satellites in the HSC-SSP than the 
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red band in Figure 9, and the discrepancy between simu- 
lations and observations could be smaller. 

Lastly, yet uncertain assumption is the isotropy of the 
satellite distribution; bright satellites are aligned in a plane 
nearly perpendicular to the Galactic plane (e.g., Kroupa 
et al. 2005; McConnachie & Irwin 2006; Pawlowski et al. 
2012; Pawlowski et al. 2013; Pawlowski et al. 2015). Hence, 
the sky-coverage correction (blue bands in Figure 9) would 
have fluctuations in the small footprint of HSC-SSP (2.896 
of the whole sky). 

In the near future, it is expected that the Vera C. 
Rubin Observatory LSST (Ivezié et al. 2019) over the en- 
tire southern sky and the Chinese Space Station Telescope 
(CSST) survey (Zhan 2011; Cao et al. 2018; Gong et al. 
2019) covering ~ 17,500 deg? will discover a significant ma- 
jority of the MW satellites plus many dwarf galaxies in and 
out of the Local Group (Qu et al. 2023). These upcoming 
surveys determine whether this new, too many satellites 


problem is real, at much higher statistical significance. 


5 Conclusions 
From the HSC-SSP Wide layer data obtained before 2021 


January (final internal data release S21A), covering ~ 
1,140 deg? in the sky, we have searched for new MW satel- 
lites using deep g, r, and i band photometry. In our previ- 
ous papers using the released data before 2018 April over 
~ 676 deg? (S18A), we reported the discovery of three can- 
didate satellites, Vir I, Cet III, and Boo IV (Homma et al. 
2016; Homma et al. 2018; Homma et al. 2019). Here in 
the data of 821A, we have re-identified these satellites and 
also found two new candidates, Sext II and Vir III, at 
Do = 126 kpc and 151 kpc, respectively. Sext II is char- 
acterized by My ~ —3.9 mag and r, ~ 154 pc, whereas 
Vir III has My ~ —2.7 mag and ra ~ 44 pc. Comparison 
with My vs. rm; relations for MW globular clusters and 
dwarf satellites indicates that these two new candidates 
most likely can be identified as UFDs. 


previously known satellites in the HSC footprint, Sextans, 


There are four 


Leo IV, Leo V, and Peg III, giving a total of nine satel- 
lites. We have presented the structural parameters of these 
satellites (except Sextans) using the current HSC-SSP pho- 
tometry data. 

To understand implications for the missing satellites 
problem in ACDM models, we have adopted the recent 
work for the expected population of MW satellites by 
Nadler et al. (2020). 


the galaxy-halo connection with simulated MW-mass host 


Their work combined a model of 


halos that associate an LMC analog and was marginalized 
to the observed satellite populations by DES and PS1 over 
^ 15,000 deg?. The expected total number of satellites 


in their model is N(My < 0) = 220+ 50 satellites with 
ra > 10 pc within the virial radius of the MW. The ap- 
plication of this model to the HSC-SSP survey after mak- 
ing a sky-coverage and completeness corrections suggests 
that a total of 3.9 -E 0.9 satellites are expected in the HSC 


footprint, compared to the nine observed satellites. This 


tension is not solved even if we exclude a classical dwarf, 
Sextans, in the galaxy count or even if we adopt the pre- 
diction of a different model for the underlying population 
of MW satellites (Newton et al. 2018). Therefore, while 
the high discovery rate of MW satellites, combined with a 
likely model of the galaxy-halo connection, appears to sug- 
gest that the missing satellites problem for ACDM models 
is no more serious than previously thought, we encounter 
a new tension, an apparent too many satellites problem. 
Meanwhile, it is important to recognize that our results are 
based on small-number statistics and several assumptions 
such as the isotropic distribution of satellites and the de- 
tection limit distance depending only on their luminosity. 

To get further insights into this issue, we require spec- 
troscopic follow-up studies of stars in these candidate 
satellites to constrain their membership as well as their 
chemo-dynamical properties. We also need to explore high- 
resolution imaging of these satellites with space telescopes 
(HST, Euclid, and more) to clearly make star/galaxy sep- 
aration for faint sources. A tighter constraint on the MW 
satellite population will be brought by the LSST survey of 
the Vera C. Rubin Observatory over the entire southern 
sky and the CSST survey covering ~ 17,500 deg”, provid- 
ing the number density of satellites and their radial dis- 
tribution at a much higher statistical significance than the 


current work. 
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Appendix 1 Space, CMDs, and density 
distributions for other seven satellites 


Here we show the space and CMD distributions for four 
known satellites, Sextans, Leo IV, Leo V, and Peg III, and 
candidate UFDs that we reported in our previous papers 
(Homma et al. 2018; Homma et al. 2019), Vir I, Cet III, 
and Boo IV (Figures 10-16). The density distributions of 


these satellites are summarized in Figure 17. 
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Fig. 10. Sextans. Same as Figure 3 except that the isochrone filter is at (m — M)o = 19.8. We note that another overdensity located at (RA, DEC)=(151.3, 
0.1)[deg] is the globular cluster, Pal 3 (Harris 1996). (c) The solid red curve shows an ellipse with a major axis of r = 0.5r;, (rj = 20.7’) and an ellipticity of 
0.25 (Tokiwa et al. 2023), whereas dotted red lines show annuli with radii r — 8.00r;, and r — 8.03r;, from the center of the galaxy. 
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Fig. 11. Leo IV. Same as Figure 3 except that the isochrone filter is at (m — M)o — 21.0. (c) The solid red curve shows an ellipse with a major axis of 
r = 2.0ra (rj = 2.6’) and an ellipticity of 0.17. 
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Fig. 12. Leo V. Same as Figure 3 except that the isochrone filter is at (m — M)o = 21.3. (c) The solid red curve shows an ellipse with a major axis of 
r = 2.0rp,(rp, = 1.2’) and an ellipticity of 0.37. 
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Fig. 14. Virgo |. Same as Figure 3 except that the isochrone filter is at (m — M)o = 19.8. (c) The solid red curve shows an ellipse with a major axis of 
r = 2.0ra (rn = 1.8’) and an ellipticity of 0.59. 
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Fig. 15. Cetus Ill. Same as Figure 3 except that the isochrone filter is at (rn — M)o = 22.0. (c) A solid red curve shows an ellipse with a major axis of 
r = 2.0ra (ri, = 1.2) and an ellipticity of 0.76. 
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Fig. 16. Bootes IV. Same as Figure 3 except that the isochrone filter is at (m — M)o — 21.6. (c) A solid red curve shows an 
r = 20r, (rj = 7.6’) and an ellipticity of 0.64. 
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Fig. 17. The density distribution of the stars in satellites passing the isochrone filter shown in Figure 10-16(d), in elliptical annuli as a function of mean 
radius. The uncertainties are based on Poisson statistics. The red line denotes the best-fit exponential profile within 5 ~ 6r;, plus a constant representing the 


background. 


